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Abstract 

We study the stability of five-dimensional Myers-Perry black holes with 
equal angular momenta which have an enlarged symmetry, U(2). Using this 
symmetry, we derive master equations for a part of metric perturbations which 
are relevant to the stability. Based on the master equations, we prove the 
stability of Myers-Perry black holes under these perturbations. Our result 
gives a strong evidence for the stability of Myers-Perry black holes with equal 
angular momenta. 



1 Introduction 



Recently, higher dimensional black holes have attracted much interest [1]. In addi- 
tion to Myers-Perry black holes [2], which have been found as the higher dimensional 
generalization of Kerr black holes, there are many exotic black holes in higher di- 
mensions [3-11]. It is quite important to study the stability of these black holes. 

The stability of higher dimensional Schwarzschild black holes has been shown [12- 
14]. For Myers-Perry black holes, the stability analysis is more difficult because of 
the difficulty of separation of gravitational perturbation equation. 1 Nevertheless, 
there have been several works on the stability analysis. In the case of Myers-Perry 
black holes with equal angular momenta in odd dimensions higher than five, it was 
shown that the special modes can be reduced to single ordinary differential equa- 
tions. For these specific modes, the stability has been shown [32]. In the case of 
five-dimensional Myers-Perry black holes with equal angular momenta, we have de- 
veloped a method of analyzing the stability by focusing on the spacetime symmetry 
£7(2) [33]. This method has been proved to be useful for the stability analysis of 
other U{2) symmetric black holes [34,35]. The purpose of this work is to study the 
stability of five-dimensional Myers-Perry black holes with equal angular momenta by 
making use of our method. We extend the previous stability analysis of Myers-Perry 
black holes with equal angular momenta [36]. 

The organization of this paper is as follows. In $2j we introduce Myers-Perry 
black holes with equal angular momenta and discuss the spacetime symmetry. In $3], 
we explain how to find master variables with which one can deduce single ordinary 
equations. In $H the master equation for most symmetric mode is derived. The 
stability for this case is analytically shown. In the master equations for other 
modes are derived and the stability of Myers-Perry black holes for these modes is 
numerically shown. The final section is devoted to conclusions. In Appendices [A] 
and [B], detailed calculations for deriving master equations are shown. 



2 Myers-Perry black holes with equal angular mo- 
menta and their symmetry 

In this section, we introduce five-dimensional Myers-Perry black holes with equal 
angular momenta and clarify the symmetry of the spacetime. 
The metric is given by 2 

ds 2 = ~dt 2 + ^ + T -^o+o- + ( ^) 2 } + 2 ± (eft + V) 2 , (1) 

1 The separability of geodesic Hamilton- Jacobi equations, Klein-Gordon equations and Dirac 
equations in the Myers-Perry spacetime have been shown in Refs. [15-31]. 

2 In our previous work [33], we used a different coordinate. To put the metric into this form, we 
need the coordinate transformation r 2 + a 2 — > r 2 and a change of the parameter m — > 2fi. 
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where 

and parameters \i and a are defined for \i > and a 2 < /i/2. Here, we define the 
invariant forms o a (a = 1, 2, 3) of SU(2) satisfying the relation do a = l/2e abc a b A o c 
as 

a 1 = — sin i/jdO + cos ■?/> sin i 
o" 2 = cos + sin ip sin I 
a 3 = dip + cos 6^0 , (3) 

and made a combination 

^ = ^T^ 2 ). (4) 

The coordinate ranges are < # < 7r, < </> < 27r, and < ip < in. The dual 
vectors of a a are given by 

e x = — sin H — cot # cos ■j/jcL , 

e2 = cosipdg + — — %ds — cot 9 smtpd^ , 

e 3 = > ( 5 ) 
and, by definition, they satisfy the relation a1e\ = 5%. 

The horizon r = r + is found by solving G(r + ) = 0, namely, 

r\ = fi + a/ fi(fi - 2d 2 ) . (6) 
The angular velocity of the horizon is given by 

n H = ± . (7) 
' + 

In term of r + and f2#, the two parameters (a,/i) in the metric ([TJ) can be rewritten 

as 

We see the the upper bound of the angular velocity, 

n H < -A- = nsr • (9) 

Apparently, the metric (pQ) has the SU{2) symmetry characterized by Killing 
vectors £ a , (a — x, y, z): 



£ x = cos 4>dg H ^cL — cot 6 sin 

S1I1& 



£„ = — sin (pdg H -cL, — cot 6> cos 

sin 

6 = ■ (10) 
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The symmetry can be explicitly shown by using the relation C^ a a a = 0, where C^ a 
is a Lie derivative along the curve generated by the vector field £ Q . From the metric 
(PQ), we can also read off the additional U(l) symmetry, which keeps part of the 
metric, o + a~ , invariant. Thus, the symmetry of five-dimensional degenerate Myers- 
Perry black hole becomes SU(2) x f/(l) ~ f/(2). 3 The Killing vectors e3, £ x , £y and 
£s constitute the symmetry, i.e., e% is a generator of U(l) and £ Q (a = x,y,z) are 
generators of 577(2). 

Let us define the two kinds of angular momentum operators 

L a = i£a , W a = ie a , (11) 

where a, /?,■■• = x, y, z and a, b, ■ ■ ■ = 1, 2, 3. They satisfy the commutation relations 

[L a , Lp] = ie^L^ , [W a , W b ] = -ie abc W c , [L a , W a ) = , (12) 

where e a p^ and e a bc are antisymmetric tensors that satisfy ei 23 = e xyz = 1. Note 
that L 2 = L 2 a = W 2 . The symmetry group, U(2) ~ SU(2) x Z7(l), is generated by 
L a and W3. Here, we should note the fact that 

Cw^^ — ic 1 * 1 , Cw 3 &^ — • (13) 

Let us construct the representation of U{2). The eigenfunctions of L 2 are degen- 
erate, but can be completely specified by eigenvalues of other operators L z and W3. 
They are called Wigner functions and are defined as 

L 2 D J KM = J(J+1)D J KM , L Z D J KM = MD J KM , W 3 D J KM = KD J KM , (14) 

where J = 0, 1/2, 1, 3/2 • • • and M, K = -J, - J + 1, • • • , J. From Eqs. (JUj), we 
see that -D^- M forms the irreducible representation of U(2) ~ SU(2) x Z7(l). The 
Wigner functions are functions of (9, <fi, ip) and satisfy the orthonormal relation: 

d9 dip dip sin 9 D J KM {0, (p, 4>)D$,* Mt (9, (p, tp) = 8jj,8 kk ,6 M m> ■ (15) 
Jo Jo Jo 

The following relations are useful for later calculations 

W+D J KM = ie K D J K _ lM , W.D J KM = -ie K+l D J K+lM , W 3 D J KM = KD J KM , 

(16) 

where we have defined W± = W 1 ±iW 2 and e K = y/(J + K)(J — K + 1). From this 
relation, we obtain the differential rule of Wigner functions as 

9+Dk M = £kD j k _ 1m , d-D J KM = -t K +iD J K+1M , d-iD J KM = -iKD J KM , (17) 

where we have defined d± = e l ±di and ^3 = e\di. 

3 The spacetime |T]) also has a time translation symmetry generated by d/dt. Because of this 
symmetry, we can separate the time dependence of fields as oc e~ lut . However, this is obvious and 
we will not pay much attention to this symmetry hereafter. 
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3 A way to find master variables 



For the stability analysis of the Myers-Perry black hole (jTJ) , it is necessary to find 
master equations for metric perturbations. Because Myers-Perry spacetime (JTJ) has 
the U(2) symmetry, a group theoretical method with a twist can be used [33,37]. 
In this section, we explain a way to find master variables. Once the master variables 
are found, we can deduce the master equations from the gravitational perturbation 
equations 

- g»u(V p V a h pa - V 2 h - R pa h pry ) - Rh pu ] = , (18) 

where V M denotes the covariant derivative with respect to g^, and we have defined 

fiW 5V + K» and h = 9 pv h pv - 

Now, we consider the mode expansion of h pv . The metric perturbation can be 
classified into three parts, h A B, h A i and hij where A,B = t,r and i,j = 9, <fi, if). They 
behave as a scalar, vector and tensor for coordinate transformations of 9, (ft and if), 
respectively. The scalar h A B can be expanded by Wigner functions immediately as 

h A B = J2 h AB(x A )D K (x*) . (19) 

K 

Here, we have omitted the indices J and M because the differential rule of Wigner 
function f fT7l) cannot shift J and M and therefore the modes with different eigen- 
values J and M are trivially decoupled in the perturbation equations. 

To expand the vector part h^i, we need an elaborate method. First, we change 
the basis {<9j} to {e a }, that is, Ha% = ^Aa^t where a = ±,3. Then, because tiAa is 
scalar, we can expand it, using Wigner functions as 

hAiix") = h A+ (x»)at + h A -{x»)ar + h A3 (x»)af 

= E [h K A+ (x A )°tD K ~i + h K A _{x A )o;D K+1 + h%(x A )afD K ] . (20) 

K 

In the expansion of h A +, Ha- and Has, we shift the index K of Wigner functions, for 
example, h A+ is expanded as J2k ^a+Dr-i- The reason is as follows. The invariant 
forms a ± and a 3 have the U(l) charge ±1 and 0, respectively (see Eq. flTSjl ). while 
Wigner function D K has the U(l) charge K (see Eq. (114]) ) ■ Therefore, by shifting 
the index K, we can assign the same £7(1) charge K to a^Dx-i, (t^Dk+i and <j\Dk 
in Eq. (l2Q|). 

The expansion of tensor part can be carried out in a similar way as 

hij(^) = Yl [ h ++°t°t D K-2 + 2h«_ojojD K + 2h* 3 <T?<T?D K -i 

K 

+h K _arajD K+2 + 2h K ? arap K+1 + h^ala]D K ] . (21) 
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To assign the same U(l) charge K to each term, we shift the index K of Wigner 
functions. 

Substituting Eqs. (TT9l) . (l20l) and (12T1) into the perturbation equations l f l8 | . we 
obtain the equations for each mode labelled by J, M and K. Because of U(2) 
symmetry, different eigenmodes cannot appear in the same equation. 

It is interesting that we can find master variables from the above information. 
First, we should note that coefficients of the expansion have different indices K and, 
therefore, coefficients of components h AB , h^ a and h^ b are restricted as follows: 



K+ 


h A +, h + s 


h-AB, 


^ A3, h + -, h 33 


h A -,h- 3 


h-. 






K-2 


< J 




K-l 


< J 




K 


< J 




K + l 


< J 




K + 2 


< J 



For example, for the J = mode, we can classify the metric perturbation as follows: 



K+ 


h<A+, h + s 


h>ABi h-As, h + -, h 33 


h A _,h_ 3 




K = 2 












K = 1 












K = 












K= -1 












K=-2 



In the above table, variables in each row can couple with each other. Apparently, 
h±± are decoupled, hence, it is straightforward to obtain the master equation for 
these variables. Other variables, (h A +, h+s), (h A B, h A3 , h + -, h 33 ) and (h A -, ^-3), are 
coupled in each set. However, after fixing the gauge degrees of freedom, we have 
the master equation for each set. In total, there are five master equations, which 
matches the physical degrees of freedom of the tensor perturbations at this level. 
We can continue this analysis. For J = 1/2 modes, we have the following table: 



K+ 


h A +, h +3 


h>ABi h A3 , h+-, h 33 


h A -, h_ 3 




K = 5/2 










# = 3/2 


K = 3/2 










K = 1/2 


K = 1/2 










K = -1/2 


K = -1/2 










K = -3/2 


K = -3/2 










K = -5/2 



From the above analysis, we see that (J,M,K = ±(J + 2)) modes are always 
decoupled. Each mode can be reduced to the single master equation. Thus, we 
have shown that we obtain an infinite number of master equations for the metric 
perturbations, although they are not everything. 

Because of the relations h*__ = h++, h* A _ = Iia+ and h*_ 3 = h +3 , we see that 
(J, M, —K) modes are complex conjugate of (J, M, K) modes. Therefore, we will 
assume K > in the following sections. 
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4 Stability analysis for ( J = 0, M = 0, K = 0) mode 



In the previous section, we showed how to find master variables. In this section, we 
will derive the master equation for the J = M = K = mode and show the stability 
for this mode. The stability for this mode has previously been shown in Ref. [36]. 
However, we will show the stability again using our formalism. 

For the ( J — 0, M — 0, K — 0) mode, we must consider metric components 
h u , htr, h rr , h t3 , h r3 , h + _ and h 33 . We set as 

hp,da?dx v = Re{e~^' [h tt (r)dt 2 + 2h tr (r)dtdr + h rr (r)dr 2 + 2h t3 (r)dta 3 

+ 2h r3 (r)dra 3 + 2h + ^(r)a + <j- + h 33 (r)a 3 a 3 ] } , (22) 

where Re represents the real part of a complex quantity. With the gauge parameters 

= Re{^(r)e-^} , = Re{£ 3 (r) e -^ 3 } , (23) 

the gauge transformations for those components are given by 

6 htt = -2zuj^ t -3 — fr , 



H 4/z 8/xa 

Ohtr = Kt- 3nl \ it ~ l^ir + £3 , 

r 6 G{r) r t) G(r) 



2fiaG(r) . 
oh t3 = 4 r - iu£ 3 . 

r t) G(r) 



s-i _ 4/xa , 2(^-2/1^-2^^ 

r d G(r) r°G(r 
<5/i + _ = rG(r% 



r 



4 



2/ia 2 )G(r) 



^33 = ^ ^T-^tr , (24) 

where ' = d r and 5 represents the gauge transformation defined by 5h^ u = + 
V^. The following gauge conditions fix the gauge degrees of freedom completely: 

h tt = 0, h t3 = 0, h 33 = 0. (25) 

Substituting Eqs. (122!) and (|25|) into Eq. (TTBl . we obtain a set of ordinary differ- 
ential equations. These equations can be found in Appendix lAl Eliminating h tr , h rr 
and h r3 from these equations, we obtain the Schrodinger-type master equation, 



dr[ 



+ V (r)$ = o; 2 $ , (26) 
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where we define the new variable 



°~ r 3 / 2 (3r 4 + 2/ia 2 ) + " ' 1 ' 



and the tortoise coordinate 



(r 4 + 2/xa 2 ) 1 /2 



The potential Vq is given by 



Vb(r) = — ^77^ [315r 20 + 162ur 18 + 2430ua 2 r 16 

ov ; 4(3r 4 + 2/ia 2 ) 2 (r 4 + 2/ia 2 ) 3 r 2 L P P 

+ 1392/i 2 a 2 r 14 + 5400/xVr 12 + 5808 / u 3 a 4 r 10 + 2608/A 6 r 8 

+ 6080/i 4 a 6 r 6 - 2064/i 4 a 8 r 4 + 32/i 5 a 8 r 2 - 160/i 5 a 10 ] . (29) 

To prove the stability, we must show the positivity of Vq. The typical profiles 
of V Q shown in Fig. CD indicate that the potential is always positive. In fact, the 
positivity can be proved from the expression (1291 . To see the positivity of Vo(r), we 
focus on r 6 , r 4 and r° terms in the big bracket in Eq. (1291 . After dividing them by 
16/i 4 a 6 , we collect them as 

/(r) = 380r 6 - 129a 2 r 4 - 10/ia 4 . (30) 

If f(r) is positive, V (r) is also positive. From Eq. (jHJ), we see // < r 2 and a 2 < r+/2. 
Therefore, 



f{r) > 380r 6 - — r 2 r 4 - |r« = 313r 6 + — r 4 (r 2 - r 2 ) + -(r 6 - r«) > . (31) 
2 2 2 2 



This proves the stability for the J = M = K = mode. 

5 Stability analysis for K ^ modes 

5.1 Master equation for (J = 0, M = 0, if = 1) mode 

As we have shown in for the ( J — 0, M = 0, K — 1) mode, we have to consider 
components h t+ ,h r+ and /i+3, namely, 

h^dx^dx" = Re{2h t+ (t,r)dta + + 2h r+ {t,r)dr a + + 2h +3 (t,r)a + a 3 } . (32) 

With the gauge parameter £j(x M ) = Re{£ + (t,r)af}, the gauge transformations for 
these components are given by 

6h t+ = i + + ^ + , 5h r+ = ? + -k + , ^+3 = ^£+, (33) 




Figure 1: Typical profiles for the potential V are depicted. From top to bottom, 
each curve represents the potential for f2#/f2™ x = 0.1, 0.7, 0.9 and 0.99. We see the 
positivity of these potentials. 



where ' = d/dt and ' = d/dr. To derive the master equation, it is convenient to use 
the action instead of equations of motion (fTSl) . The action for metric perturbation 
in a vacuum is given by 

+ 2V A P V^ W " 2VVV"/i] , (34) 

where we use the unit 167rG 5 = 1; here G5 is the five-dimensional Newton's constant. 
Now, we define gauge invariant variables, 

„2 



h t . 



ir 



ft = — 5- + ~ n^+3 n^+3 5 fr = — T- * 

ar z 



h r 



r+ 



2fia* ar* 
In term of these variables, the action f|34l) becomes 



2/ia 



(35) 



S 



dtdr 



r(r 4 + 2/ia 2 ) • 2 4/i 2 a 4 2 

— I J r- — Jt I + K^y/ \ \ Jt\ 
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r 5 G(r) 



+ V ° 2(r 5 2 fl2) |M 2 + 2/mIm{r(/ r - / t ')/ r * - 4/J r *} 



• (36) 



In the above action, there are two fields f t and / r . However, f t is not the physical 
degree of freedom. Therefore, we can eliminate it from the action and get the master 
equation. In Appendix \tQ we show the details. As a result, in terms of the new 
variable 



$1 



„3/2 / ^,4 



2/ia 



2\l/4 



(r 10 + 2/xa 2 r 6 + /i 2 a 6 )V2 



r(r 4 + 2/ia 2 



■(/r - /t) +i^arf r 



(37) 
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the master equation can be obtained as 



drl 



+ y 1 (r)$ 1 = [^-2n 1 (r)] 2 $ 1 , (38) 



where 



n ^ _, 2 /* a A « 2 ^ 4 (5r 4 + 6/xa 2 )G(r) 
U; r 4 + 2/ia 2 V 4(r 10 + 2/ia 2 r 6 + /i 2 a 6N ' ' 1 



and 

VUr) =^r— A ONQ/ G ^ — [35r 32 + 18ur 30 + 310ua 2 r 28 

1{ - J 4r 2 (r 4 + 2/ia 2 ) 3 (r 10 + 2/ia 2 r 6 + / u 2 a 6 ) 2L P P 

+ 160^ 2 a 2 r 26 + 1192/xVr 24 + 2 / u 2 a 4 (152/i - 75a 2 )r 22 + 3068/i 3 a 6 r 20 

- 64/i 3 a 6 (2/i + 15a 2 )r 18 + 5208//a 8 r 16 - 16/i 4 a 8 (30/i + 133a 2 )r 14 
+ 3/i 4 a 10 (1424/i + 5a 2 )r 12 - 1654/iW + 2 / uV 2 (432 / u + 25a 2 )r 8 

- 168/i 6 a 14 r 6 + 68/iW - 24/iVV + 56//a 18 ] . (40) 

We used the tortoise coordinate defined in Eq. (l28l) . In Eq. (l38l . we separated the 
t dependence of $i by the Fourier transformation $i(t, r) = e~ lwt <&i(r). 
The asymptotic forms of fii(r) and Vi(r) are 

fil ( r ) (r -> oo) , fix(r) -> fi H , (r -> r + ) (41) 

and 

Vi(r)->0, (r^r+,oo) (42) 

where is the angular velocity of the horizon defined in Eq. (|7j). Therefore, the 
asymptotic form of the solution of master equation (1381 becomes 



$! -> e ±4wr * , (r -> oo) , $i -> e ±4{w " 2Qff)r * . (r -> r + ) (43) 

Before discussing the stability for this mode, we shall derive the master equation 
for (J, M, K = J + 2) modes. The stability analyses of both cases will be discussed 
simultaneously. 



5.2 Master equations for (J, M, K = J + 2) modes 

To obtain the master equations for (J, M, K = J + 2) modes, we set /i^ as 

h^x^dx^dx" = Re[h ++ (t,r)D J (x i )a + a + ] , (44) 

where L>j = D J K=J M . This field is gauge invariant. Substituting Eq. (1441) into 
Eq. (l34l) and using the differential rule of Wigner functions (fT7j ). we obtain the action 
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for the K = J + 2 mode, 



S = — I dtdr 
4 



r 4 + 2/m 2 
4r 5 G(r) 
1 



1/^ 



G(r) 
4r 



i(J + 2)ixa 
r 5 G(r) 



{h ++ h\ 



{ - (J + 1)(J + 2)r 8 + 2/i( J 2 + 3 J + l)r 
+ 2/x(2/x + ( J + 4)a 2 )r 4 - 12/iVr 2 + 8/iV} |/i 



++I 



(45) 



where the asterisk denotes the complex conjugate. We can derive the equations 
of motion for (J,M,K = J + 2) modes from the above action. Defining the new 
variable 

^+2 = ' J/2 h ++ > ( 46 ) 



we obtain the master equation, 



d 2 $ 



J+2 



drl 



+ V J+2 {r)Qj + 2 = [to- 2( J + 2)fi J+2 (r)] 2 $ J+2 



where 

and 

Vj +2 (r) 



2fia 



r 4 + 2/ia 2 



(47) 



(48) 



G(r) 



, rl(r , 2/ , f|2) 3 K 4J + 7 )( 4 J + 5 ) r + 18 ^ r + W(16 J 2 + 32 J + 5)r s 
- 40/iVr 6 - 4 / u 2 a 4 (16J + 35)r 4 + 8/?a 4 r 2 - 40/i 3 a 6 ] . (49) 

We separated the t dependence of $j+2 by the Fourier transformation $j +2 (t, r) = 
e-^$ J+2 (r). 

Since the asymptotic forms of nj +2 (r) and Vj +2 (r) become 



JWr)->0 (r^oo) 



^j +2 (r) 



Qjj , (r -> 



and 



V> +2 (r)->0, (r^r + ,oo) 
we get the asymptotic form of solution of master equation fT4"7l) as 



J+2 



-,+iwr. 



r — > oo 



J+2 



\r — ► rj 



(50) 
(51) 

(52) 
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5.3 Method to study the stability of K ^ modes 

Since the master equations for K ^ modes are not of the Schrodinger type, we 
cannot show the stability for these modes from the positivity of the potential as was 
done in the case of the J = M = K = mode. Here, we will follow the method 
used to show the stability of Kerr black holes [32,38]. 

To discuss the master equations for (J = 0, M = 0, K = 1) and (J, M,K = J+2) 
simultaneously, we write the master equations as 

" d ~TY- + V K (r)$ K = [u- 2KQ K {r)f<$> K . (53) 

For K — 1 and K = J + 2, the above master equations reduce to Eq. (1381) and 
Eq. (1471) . respectively. 

Now, we will examine if there exists a quasi-normal mode with Imw > in 
the system. Because of the time dependence h^ v oc e~ luJt , the existence of such a 
mode implies instability of the system. Recall that the boundary conditions for 
quasi-normal modes are given by 

®k -> e - i{uj - 2KnH)r * , (r -> r+) 

$ K - Z oat e^ r * . (r -> oc) (54) 

Namely, the wave function must be ingoing at the horizon and outgoing at infinity. 

To study the stability using the master equation f[53l) , we start with the assump- 
tion that Myers-Perry black holes are stable for sufficiently small angular velocity 
Qh- This is a natural assumption because the higher dimensional Schwarzschild 
black hole is stable [13]. In fact, in the limit of a — > 0, the master equation (ioBl 
takes the Schrodinger form and the positivity of the potential is easily seen. Under 
this assumption, for small Qjj, the imaginary part of the quasi-normal frequency 
must be negative, Imw < 0. Now, if there exists instability, a quasi-normal mode 
with Imw > will appear at some point as we increase Qh- This means that one 
of quasi-normal modes must cross the real axis in the complex to plane for some 
Qh- Therefore, if the black hole is unstable for large Qh, there must be a critical 
value Qh = H™' for which there exists a mode with Imw = under the boundary 
condition (|54l) . We look for such Q ^. 

For the numerical analysis, it is convenient to define quasi-normal modes as 
follows. Given the ingoing wave at the horizon, 

we generally obtain the wave function at infinity as 

$x - Z ont e^ + Z iQ e-^ r * , (r - oo) (56) 

where Z out and Z in are constants. Hence, quasi-normal modes are defined by the 
condition Z- m = 0. 
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For the purpose of searching for fl™ 1 , we take u to be real. In this case, the 
Wronskian of &k is conserved, that is, 



Im 



, (57) 



for any r\ and r 2 . We take T\ = r + and r 2 = oo. Then, from Eq. (1571) . we obtain 
the relation, 

2K{l H - lu = cu\Z out \ 2 , (58) 

where we used the boundary condition (1541 ). Since the master equation ( 1531) is 
invariant under uo — ► — and a — > —a, we take > 0. Then, from Eq. (l58l) . we 
obtain the inequality 

< w < 2Xn H . (59) 

We should recall that VLh must also satisfy Eq. ((9)). We will search for Q^ 11 in this 
region. 

5.4 WKB analysis 

Before the numerical analysis, we study the master equation (153]) using the WKB 
approximation. Let us define 

V K {r) = V K (r) - (u - 2Kn K (r)) 2 . (60) 

If Vk < everywhere, one of the WKB solutions of the master equation ( |53l) is 



®K ~ exp \ i J dr,\/-V K ) . (61) 

The asymptotic form of this solution becomes 

®K -> e - i( "- 2 ™ ff)r * , (r -> r+) 

$ K -> e iu}r * , (r -> oo) (62) 



where we use Eq. (1591 . Equation (1621 ) nothing but the boundary condition (I54I) . On 
the other hand, if there is a region satisfying Vk > 0, the WKB analysis leads to 
\Z 0Ut /Z- m \ ~ 1. Thus, the condition 

V(r) < for all r (63) 

can be considered as a rough criterion for the existence of instability. Let us look at 
the potentials V\ and V 2 shown in Fig. [2l The potential V 2 tends to be negative for 
large Qh and u. For K > 2, we find a similar behavior to that of V 2 . In the case of 
Vi, however, a positive region remains even for sufficiently large u> and Qh- From 
the results of these WKB analyses, we can speculate that the (J = 0, M = 0, K = 1) 
mode is stable and the (J, M, K = J + 2) modes might be unstable for large Qh and 
io. With this intuition, we shall perform numerical calculations. 
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1 1 5 2 2.5 3 3.5 4 4.5 5 1 1.5 2 2.5 3 3.5 4 4.5 5 

r/r + r / r + 

Figure 2: V\ and V 2 for f2#/f2g ax = 0.99. We see that V 2 tends to be negative for 
large u. However, for Vi, a positive region remains even for sufficiently large u. 

5.5 Numerical analysis 

Now, we must solve Eq. (l53l) numerically. We integrate this master equation from 
ri = (1.0 + 10~ 6 )r + to r 2 = 100. 0r + with the initial condition <3>x = exp(— iujr*) 
at r = r\. This is nothing but the ingoing boundary condition at the horizon. At 
r = r 2 , we checked the ratio of the amplitudes of ingoing and outgoing modes: 

Z(u,n H ) = \Z out \ 2 /\Z in f . (64) 

We calculated this ratio for each u and fin in domains j9]) and ( 1591) . If Z in = at 
r = r 2 for some u and the function Z(u, fin) would diverge. This would be a 
signal of instability. We plot 

Z max (tt H ) = max Z(u, Q H ) (65) 

fixed Uh 

in Fig. [3l For (J = 0, M = 0, K = 1), we see Z max ~ 1. It is expected from the 
WKB analysis in the §5.41 For (J, M, K = J + 2), the figure shows that Z max takes 
a maximal value at Qh = fi^ ax - This behavior is also expected from the results of 
the WKB analysis in the previous subsection. Taking a look at Fig. [3l at least for 
Q H < f^g ax , Z max (Q H ) is finite. 4 Therefore, we conclude that Myers-Perry black 
holes are stable when Qh < fi^ ax . 

6 Conclusions 

We have studied the stability of five-dimensional Myers-Perry black holes with equal 
angular momenta. Utilizing the symmetry of these solutions, SU{2) x U(l) ~ U(2), 

4 Note that, in the maximally rotating case, the outer and inner horizons are degenerate and 
the spacetime structure changes from that in the fin < fl^ ax case. Therefore, we must analyze 
the case of fin = fl^ ax separately. 
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Figure 3: Function Z max (n H ) for ( J = 0, M = 0, K = 1) mode and (J, M,K = J+2) 
modes. We see that Z max is finite for Qh < ^# ax . 



we have identified the master variables with which the perturbed equations can be 
reduced to master equations. By analyzing the master equations, we have shown the 
stability of Myers-Perry black holes with equal angular momenta. Strictly speaking, 
we have not shown the stability of Myers-Perry black holes completely, because we 
analyzed restricted modes. Empirically, however, the instability appears in the lower 
eigenvalue modes. For example, the Gregory-Laflamme instability appears in the 
s-wave. Therefore, our results for ( J = 0, M = 0, K = 0, ±1, ±2) modes give strong 
evidence for the stability of Myers-Perry black holes. 

In this paper, we considered the Myers-Perry black holes with equal angular 
momenta. There are several qualitative arguments for the stability of general Myers- 
Perry black holes. In the case of D > 6, the argument goes as follows. The Myers- 
Perry black hole with one rotating axis reduces to a black brane solution in the 
limit of large angular momentum. Since Gregory-Laflamme instability occurs in the 
black brane system [39], we can speculate that the Myers-Perry black hole with one 
large angular momentum is unstable in the (D > 6)-dimension [40]. In the case 
of D = 5, we cannot use the above argument. Instead, we can use the fact that 
the black ring solution is thermodynamically more preferable than the Myers-Perry 
black hole with one large angular momentum in order to infer the instability of 
the five-dimensional Myers-Perry black hole with one large angular momentum [4]. 
Although these arguments seem to be reasonable, the direct dynamical analysis must 
be performed to achieve the final conclusion. 

The Myers-Perry black holes can be generalized to those with the cosmological 
constant, the so-called Kerr-AdS spacetimes [41-43]. Our method is also applicable 
to five-dimensional Kerr-AdS black holes with equal angular momenta. It is interest- 
ing to investigate such a case from the point of view of the AdS/CFT correspondence. 
From the thermodynamics of CFT4, which is a dual theory of five-dimensional Kerr- 
AdS spacetime, instability can be expected for Kerr-AdS spacetime [43,44]. On the 
gravity side, there are several works on the stability of Kerr-AdS spacetime. In 
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the case of 4-dimensional Kerr-AdS spacetime, the superradiant instability has been 
found [45,46]. In odd dimensions higher than five, the same instability of Kerr-AdS 
black holes with equal angular momenta has been shown to exist [32]. In the case of 
a D > 7-dimensional Kerr-AdS black hole with one rotating axis, it has been shown 
that the superradiant instability appears in the tensor-type perturbation [47]. How- 
ever, there is no work on the stability analysis of five-dimensional Kerr-AdS black 
holes (except for a massless case [48]). It is of interest to study the stability of 
five-dimensional Kerr-AdS black holes using our formalism. 
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A Gravitational perturbation equations for J = M = 
K = mode 

In the following, we list the components of the Einstein tensor for the ( J = 0, M = 
0,K = 0) mode. 



_ 4icun 2 a 2 , 8iufi 2 a 2 (r 2 - /i) 

o^tt — s h t n tr 

r o r a G(r) 



G(r)(3r 6 - 12«r 4 + 2/i(6« + a 2 )r 2 - 8/i 2 a 2 ) 

+ ^7 V 

+ [3r 10 - 6r 8 /i + 2/i(-6/i + fia 2 u 2 - 4a 2 )r 6 

+ 8/i 2 (3« + 7a 2 )r 4 - 4a 2 « 2 (22/i + 3a 2 )r 2 + 56a V 3 )] Kr 
4icu/ia(r 2 - 2/i) 4iu;/ia(r 4 - 4/ir 2 + 2/i(2/i - a 2 )) 

2G(r) (r 2 - 2/i) , _ 4/i(r 4 - 2(/i + 2a 2 )r 2 + 6/ia 2 ) , 

^9 h +- 

. 8/ia 2 ((/ico 2 - 2)r 6 + 2/ir 4 + 4/i(/i - a 2 )r 2 - 4a 2 /i 2 ) 

+ ' (66) 

iu;(3r 4 — 6«r 2 + 2/ia 2 ) 2u 2 fia 
dGtr = ^5 hrr ~ r *G(r) 

1^- ""^ , (67) 

_ -iufia(r A + 2/ia 2 ) , 2iufj,a(r 4 + 2/ia 2 ) (r 2 - /i) 
dG * 3 - ^6 ^ r 9 G(r) ^ 

//aG(r)(2r 4 - 3/ir 2 + 2/ia 2 ) /ia f 2 10 8 

+ 2/i(aV - 6)r 6 + 8/i(3« + 5a 2 )r 4 - 88/i 2 a 2 r 2 + 56« 2 a 4 )]/i rr 
2^(r 6 -2/ir 4 + 2/ia 2 r 2 -8/i 2 a 2 ) ^ 8 6 

^6 ^ 3+ 2G^V [ dr +i2/ir 

- 4/i(3« + a 2 )r 4 + 24« 2 a 2 r 2 - 16// 3 a 2 + 4/i 2 a 4 ]/i r3 
2fxaG(r) 2/ia(r 4 + 2/ir 2 - 6/xa 2 ) , 

+ "TZ^rT^ 10 - 6r 8 + 2/x(aV + 6)r 6 
r i4 G(r) 

- 16«a 2 r 4 + 8« 2 a 2 r 2 - 8/i 2 a 4 ]a+_ , (68) 
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_ ioo(3r 4 + 2/ja 2 ) (3r 4 - 2/ta 2 ) Aiuo^a 
r a G(rJ r° r°G(r j ) 

+ iJ $GW"'- + 7^ [ " 2r '° " 2r " + 2M<l2 " 2 + 4)r6 

- 8/i(// + a 2 )r 4 + 16/tVr 2 - 8/i 2 a 4 ] , (69) 

_iuna u; 2 (r 4 + 2/ia 2 ) 4iu/ia 
6Gr3 ~— Kr ~ 2r 4 G(r) " P^) ' (?0) 



za^r* + 2/ia 2 ) , itu(r 4 — /tr 2 + 2/ta 2 )(r 4 — 2/ta 2 ) 



- G(r)( o r2 ^ r + A^ V ° - 2 ^ 8 + 2 M«V - 2)r 6 
2r 4r 8 

+ 8/i(/i + a 2 )r 4 - 8/i 2 aV - 8/i 2 a 4 ]a 

, 2ja^a 4^/i 2 a(r 2 - 2a 2 ) 1 

+ — K 3 + Tn~T\ hr3 + o G ( r ) h +- 

r z r'G(r) 2 

(r 4 - 6/ir 2 + 10/ia 2 ) 1 r 2 10 , „ 8 

H in . [a; r + 4r 



rr 



2r 5 +- 2r 10 G(r) 

+ 2/t(aV - 8)r 6 + 16/t(/t + 2a 2 )r 4 - 64u 2 a 2 r 2 + 48/iV]a+_ , (71) 



iu{r 4 + 2/ta 2 ) 2 icj(r 2 - /i)(r 4 + 2/ta 2 ) 2 
G(r)(r 8 - /tr 6 + 4/ia 2 r 4 - 6/t 2 a 2 r 2 + 4/t 2 a 4 ) , 

+ _L LA 14 - 2r 12 + 4//(aV - l)r 10 

+ 4/t(2/i + 5a 2 )r 8 + 4/xV(aV - 16)r 6 

+ 8//V(6/i + 13a 2 )r 4 - 176/t 3 aV + 112/t 3 a 6 ]a rr 
iiOjj,a(r 4 + 2/ta 2 ) , 2ioojja(2r 6 — 3/tr 4 + 4/ta 2 r 2 — 2/t 2 a 2 ) 

^ ^ 3 + ^G(0 r3 



G(r) (r 4 + 2/ta 2 ) r 8 - 6/xr 6 + 8/mV - 4/z 2 a 2 r 2 + 12/x 2 a 4 , 

+ 2? + - 2? + - 

+ o TI777 \ [^ 2 ^ 14 - 2r ' 2 + MaV - l)r 10 + 4/t(4/t - 3a 2 )r 8 
2r i4 G(rJ 

+ 4^/i 2 a 4 r 6 - 24/xVr 4 + 16/AV - 16/t 3 a 6 ] a+_ . (72) 

From 5G tr = 6G rr = 5G r s = <5G + _ = 0, We can eliminate h tr ,h rr and h r3 , and 
obtain the master equation ( |2"6"1 ). 
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B Derivation of master equation for ( J = 0, M = 
0, K = 1) mode 

First, we should note there is a nondynamical variable f t in the action (136ft . This 
means there are constraints. To treat the constraints properly, it is convenient 
to adopt a Hamiltonian formalism [49]. From the action (l36l) . one can read off a 
Lagrangian as 

C=e A \f r -f t \ 2 + e B \f t \ 2 + e c \f r \ 2 

- *e D {(f r - fl)f: - (f r - flTfr} + ie E {f t f* r - f t *f r ) . (73) 

Here, to simplify the calculation, we have introduced the notations 

r(r 4 + 2/ia 2 ) B 4/j, 2 a 4 



4 ' r 5 G(r 
c _ V«V 2 - ^) n „ k 



fiar , e = 4/xa . (74) 



The conjugate momenta of f r , f* are given by 

*r= 7 -r = e A Ur-f t )+ie D fr, 
df* 

K =|x = Afr - flY - ie D f? ■ (75) 

df r 

The above momenta satisfy the following canonical commutation relations, 

{f r (t,r)X(trf} = {f:(t,r),7r r (t,r')} = 5(r - r') , (76) 

where { , } represents a Poisson bracket. Taking the variation of the Lagrangian 
(|73l) with respect to /*, we obtain the constraint 



ft = e- B {-n' r + ie E f r ) . (77) 
The Legendre transformation gives the Hamiltonian 

'Ho =n r fr + Kfr - £ 

=e"V< - e B f t fl - (e c - e 2D ~ A )f r f: + ir r f? + Kf't 

+ ie D - A {vrf; - Kfr) - ie E (f t f: - fifr) ■ (78) 



Using the constraint equations, we obtain the physical Hamiltonian 

^phys 



n phys = e"V< + e- B 7rX* + (e 2E ~ B - e c + e 2D ~ A )f r f: 



+ ie E - B Xf* r - TT^/r) + ^°~ A Xf: - Kfr) 



(79) 
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where total derivative terms are omitted. We obtain the following equation of motion 
for the physical variable: 

f r = e- A 7r r -(e- B Tr' r y + i(e E - B f r y-ie D - A f r , 

^ = -( e ™-B_ e C + e 2D-A )fr _ ie E-B 7r , r _ ie D-A 7rr {gQ) 

We obtain the master equation for the ( J = 0, M = 0, K = 1) mode by eliminating 
the variable f r from Eq. (1801 ). 
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